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The spectrum of photons arising from WIMP annihilation carries a detailed imprint of the struc- 
ture of the dark sector. In particular, loop-level annihilations into a photon and another boson can 
in principle lead to a series of lines (a WIMP forest) at energies up to the WIMP mass. A specific 
model which illustrates this feature nicely is a theory of two universal extra dimensions compacti- 
fied on a chiral square. Aside from the continuum emission, which is a generic prediction of most 
dark matter candidates, we find a "forest" of prominent annihilation lines that, after convolution 
with the angular resolution of current experiments, leads to a distinctive (2-bump plus continuum) 
spectrum, which may be visible in the near future with the Fermi Gamma-Ray Space Telescope 
(formerly known as GLAST). 



I. INTRODUCTION 

Cosmological and astrophysical observations 
provide compelling evidence for dark matter 
(DM), but so far they have not provided a 
smoking-gun indication of its identity. The hunt 
for dark matter is now entering a new era, with 
the current generation of indirect and direct detec- 
tion experiments closing in on weak scale masses 
and couplings, and the LHC turn-on just around 
the corner. Among the plethora of dark mat- 
ter candidates. Weakly Interacting Massive Parti- 
cles (WIMPs) appear particularly appealing, since 
they combine the virtues of weak scale masses and 
couplings, stability, and connection to the mystery 
of electroweak symmetry-breaking ! l] . 

If WIMPs do indeed have weak scale masses and 
interactions, they are expected to produce observ- 
able effects in colliders and astroparticle experi- 
ments. In particular, indirect detection experi- 
ments seek secondary particles produced by the 
annihilation of WIMPs in regions of our galaxy 
where the DM density is high, e.g. the Galactic 
center, or DM clumps. Among secondary parti- 
cles, gamma-rays play a special role, since they 
have the advantage of traveling in straight lines 
and without sizable energy losses in the local uni- 
verse, thus tracing the distribution of DM in a 
straightforward way. 

The annihilation of WIMPs into photons typi- 
cally proceeds via a complicated set of processes, 
and the final spectrum actually contains a detailed 
imprint of WIMP annihilation that can in princi- 
ple reveal features such as the WIMP spin and/or 
other particles in the dark sector. Tree-level anni- 



hilation of WIMPs into quarks or leptons (or heav- 
ier states which decay into them) provides a con- 
tinuum of photon energies, with an upper cutoff at 
approximately the WIMP mass. This continuum 
emission is in general rather featureless, but some 
models do exhibit spectacular features, such as a 
high-energy rise due to final state radiation 0. 
Furthermore, loop-level annihilation into a photon 
and X results in a line at energy 

where moM is the WIMP mass and Mx the mass 
of the other boson in the final state (in super- 
symmetric theories, X is either another photon 
or a. Z boson [^). Line emission typically has 
smaller magnitude than continuum emission, but 
the line provides a feature that helps discriminate 
against backgrounds. More generally, it may also 
be that the tree-level annihilation is into photon- 
unfriendly modes, and the line(s) may turn out to 
be prominent features in the photon energy spec- 
trum. 

At energies in the 10 GeV - 10 TeV range, 
the energy resolution of current gamma-ray tele- 
scopes such as the Fermi LAT or Air Cherenkov 
Telescopes like CANGAROO, HESS, MAGIC and 
VERITAS, is of order AE/E - 0.1. When the 
mass of X is small compared to the WIMP mass, 
different lines may merge together because of finite 
detector resolutions. When X has a mass which is 
appreciable compared to the WIMP mass, the lines 
can be distinct and separately measurable. When 
distinct, they provide information about physics 
beyond the Standard Model. For example, the 
position of the 77 line measures the WIMP mass. 
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The relative strengths of the 77 and 7Z hnes are 
a measure of the WIMP couphng to SU{2) sin- 
glets versus doublets. An unsuppressed line 
(where H is a Higgs boson) would be an indica- 
tion that the WIMP is not a Majorana fermion or 
a scalar. Finally, if there are particles in the dark 
sector whose mass is less than 2mDM, then a pair 
of WIMPs may be able to annihilate into one of 
these exotic heavy states and a photon. The ex- 
istence of such a line could be the first indication 
that such a heavy partner is present. Taken all 
together, there could be a "forest" of lines associ- 
ated with WIMP annihilations, and they represent 
a wealth of information about the theory of dark 
matter. 

A specific model which illustrates these fea- 
tures nicely is the "Chiral Square" [1], a six di- 
mensional model with universal extra dimensions 
(UED) 3 - extra dimensions in which the en- 
the Standard Model (SM) propagates. UED the- 
ories contain WIMP candidates because of rem- 
nants of the extra-dimensional spacetime symme- 
tries which forbid the lightest of the Kaluza-Klein 
excitations of the SM from decaying. The Chiral 
Square is an intrinsically six dimensional construc- 
tion, whose dark matter candidate, the "spinless 
photon" Q (the KK mode of the hypercharge vec- 
tor boson whose spin is pointing in the extra di- 
mension) has distinct phenomenology compared to 
the five dimensional case [3, H, S Il3| . Six dimen- 
sional implementations of the SM have a mech- 
anism which automatically suppresses dangerous 
higher dimensional operators leading to proton de- 
cay 111 and motivate the existence of three gener- 



jacent sides of the square are identified with each 
other, 



ations of fermions through gauge anomaly cancel- 
lation [l3|. 

The article is laid out as follows. In Section [Til 
we (briefly) review the chiral square model, and 
what is currently known about its dark matter 
candidate. In Section lllli we calculate the con- 
tinuum gamma ray spectrum from WIMP anni- 
hilations in the galaxy. Section IIVI outlines the 
calculation of the cross sections for the annihila- 
tion processes of two WIMPs into states with one 
or two final primary photons. In Section |V] we 
discuss the prospects for detection of the astro- 
physical signal and the related uncertainties. We 
finally reserve Section IVTl for the conclusions. 



II. THE CHIRAL SQUARE 

The chiral square is a model of two compact 
universal extra dimensions. The extra dimensional 
coordinates may be represented by a pair of points 
(x^, x^) living in a square region with sides L. Ad- 



(2/,0) = (0,y) 



(2/,L)^(L,y) (2) 



which is equivalent to taking a square, folding it 
along a diagonal, and smoothly gluing the edges 
together. The folding leaves the two corners of 
the square which lie along the fold (at (0, 0) and 
{L,L) ) invariant, and identifies the remaining two 
corners (at (0, L) and (i,0) ) as the same point. 

The Kaluza Klein (KK) modes of SM fields are 
labelled by a pair of integers, (j, k) which satisfy. 



> , j > 1 - 4,0 • 
Scalar fields have wave functions. 



(3) 



(4) 



C± 



jx^ ± kx^ 



with masses given (up to boundary term effects 
described below) by. 



Mn 



(5) 



where Mq is the mass of the "zero-mode" field, 
which carries no momentum in the extra dimen- 
sional directions, and has wave function /("-O) = 
1/{2L). We identify the zero modes with the Stan- 
dard Model fields. 

The KK modes of fermions are Dirac particles. 
The (4D) chirality associated with the zero mode 
works essentially as outlined above, but the oppo- 
site chirality pieces have wave functions which are 
phase-shifted. The full details are presented in [3| , 
and are not essential for our purposes. The gauge 
bosons decompose into 4D vectors and two 
4D scalars, and V^. One linear combination of 
and is eaten, level by level, by the vector 
KK modes to provide their longitudinal degrees of 
freedom. The other linear combination of and 
are physical gauge adjoint scalars in the 4D ef- 
fective theory. We follow the usual convention and 

denote them by V^''^'' . 

The residual space-time symmetry causes the 
lightest (1,0) KK particle (LKP) to be stable. 
The allowed "large" interactions can be roughly 
understood from higher dimensional momentum 
conservation, but with the important observation 
that since moving across a boundary rotates the 
direction of the momentum by 90°, the momen- 
tum in the x^ and directions are not actually 
distinct from one another. The full momentum 
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conservation, even in this folded sense, is broken 
by terms which Hve on the corners of the square 
[13| (and whose presence is required when the the- 
ory is renormahzed p3 , which nevertheless 
preserve a bulk symmetry (invariance under re- 
flections through the middle of the square), pro- 
vided the UV physics is such that equal bound- 
ary terms live at (0,0) and {L,L). Even in the 
presence of such terms, there is a Z2 KK-parity 
U under which states that have j + k equal to 
an even integer are even, and states with 7 -|- k 
odd are odd. Just as in the 5D case [H, [la| the 
boundary terms have the effect of breaking degen- 
eracies within a given KK level and introducing 
KK-number violating (but KK-parity preserving) 
couplings among the KK modes p^ . 

We follow the usual reasoning, and assume 
that the boundary terms are dominated by ra- 
diative contributions from the bulk physics [Tsj . 
While not strictly required by any theoretical or 
phcnomenological argument,^ loop- level boundary 
terms motivate having the colored and charged KK 
modes being heavier than their neutral counter- 
parts (which, at least at the (1,0) level is essential 
for the theory to contain a viable WIMP). 

Motivated by the expectations of radiative 
boundary terms, the hghtest (1,0) mode is ex- 
pected to be the scalar partner of the hypercharge 
boson, 5^^'°^ (which we will often abbreviate as 
simply Bh) ^M- While electroweak symmetry- 
breaking mixes i?^'"'' with its SU{2) counterpart, 
T^i/'°\the mixing angle is typically small, and the 
LKP is, to good approximation, pure Bh- This 
implies that the LKP is a gauge singlet, whose 
coupling to SM fermions and Higgs is controlled by 
the U{1) gauge coupling gi and the hypercharge of 
the matter field. As a further consequence of small 
boundary terms, its largest couplings are expected 
to be those which are KK number conserving (in 
the folded sense above). 

The LKP, as a real scalar, has a suppressed anni- 
hilation rate to light fermions in the present epoch 
and therefore annihilates predominantly into pairs 
of electroweak bosons WW and ZZ, Higgs bosons 
HH, and (if heavy enough) pairs of top quarks 

[6|. As a result, its thermal relic density is 
very sensitive to its mass and the mass of the zero 
mode Higgs, which provides a resonant annihila- 
tion channel when Mh ^ 2™^^ . The relic density 



^ However, analysis of large boundary terms in 5D does 
reveal that when they are large, it is generically more 
difficult to fit precision EW data fT7| . 
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FIG. 1: Continuum photon spectra, dN^/dx for two 
mass choices of Bh- 



is consistent for masses in the range of roughly 200- 
500 GeV, with narrow windows of Higgs masses 
appropriate for each LKP mass. We use the relic 
density as a rough guideline to the most interesting 
range of LKP masses for indirect detection, but do 
not strictly assume it applies; it could be that the 
LKP is not a thermal relic, or our extrapolation 
of cosmology to early times is flawed by imperfect 
understanding. 



III. CONTINUUM 7-RAY EMISSION 

As mentioned above, pairs of LKPs annihilate 
predominantly into pairs of electroweak bosons 
WW and ZZ, Higgs bosons HH, and (if heavy 
enough) pairs of top quarks tt. Far above the tt 
and HH thresholds, the annihilation fractions are 
roughly 50% BhBh ^ WW, 25% BhB_h ZZ, 
25% BhBh HH, with BhBh tt subdom- 
inant. Continuum gamma ray emission is largely 
the result of radiation from charged leptons which 
result when these massive objects decay, or decays 
of TT^s from the hadronization of strongly inter- 
acting decay products. The result is a rather soft 
spectrum of gamma rays, reminiscent of neutrali- 
nos in supersymmetric models. 

We compute the continuum spectrum using the 
micrOMEGAs code [ll|, based on a CalcHEP 
[l9| (partial) implementation of the chiral square 
model [l^l . From Fig. [l] where the differential flux 
(per steradian) from the continuum spectrum is 
plotted for Mb„ = 250 GeV and for Mb« = 450 
GeV it is interesting to point out that the spec- 
trum sharply decreases well before the value of 
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Mbh- This is a distinctive feature of this spe- 
cific chiral square model and is related to the fact 
that the dominant annihilation channels are into 
photon unfriendly modes, consisting of massive 
(and often neutral) particles which are unlikely 
to radiate high energy photons. This leaves the 
bulk of the photons coming as radiation (or af- 
ter hadronization, decays of tt^s) from the even 
softer decay products of the particles produced in 
the primary annihilation, which have less energy 
available for the final state photons. 



IV. CROSS SECTIONS FOR 7-RAY LINES 

Let us now consider the direct annihilation of 
LKP pairs into photons. Specifically, we will con- 
sider the process BhBh — > where V can be ei- 
ther a photon, Z boson or a KK excitation of a SM 
gauge field. Since the LKP has no direct coupling 
to the photon, the leading-order contributions to 
these processes occur at one loop. Examples of 
Feynman diagrams are shown in Fig. [S] Here we 
assume that the dominant contributions are from 
loops of SM charged fermions (£) and their corre- 
sponding pairs of KK partners (Cl^])- Note that, 
at a given KK level, the fermionic field content is 
doubled as compared to the SM such that each SM 

fermion has both a sing let (^i^^) and doublet (C^^^) 
KK partner. 

The matrix element for the generic process 
Bh{pi)Bh{p2) -ff'ipA)V''ipB) takes the form: 

e^/{pAyB*{pB)M^''{puP2,PA,PB) , (6) 



are the polarization tensors of 



where and 
the photon and V gauge boson, respectively. 

In general, the amplitude Ai'^'^ can be expressed 
as a linear combination of tensor structures built 
from the external momenta and the metric tensor 



g'"'. Considering the transverality of the polariza- 
tion tensors, ga'Pa = ^B'Pb — 0, the most general 
tensor structure is given by: 

M^" = Ai + B, p'^p'i + B2 pl^p^2 + B3 ptp'i 

+ Bi pIp^^ + i?5 p'Xp'i, + i?6 ptPA + Bj p'ip'i, 

+ Bs P2PA + Bg P2PB ■ 

However, WIMPs are assumed to be highly non- 
relativistic (with typical velocities of w ~ 10~^), 
such that the incoming momenta of the i?//'s 
are well approximated hy pi = p2 = p where 
P = (wb^,0). Thus, using conservation of mo- 
mentum {2p>^ = pj^ we can eliminate many 
of the above terms. In addition, since the WIMP's 
are annihilating nearly at rest, the final state prod- 
ucts are emitted back-to-back such that (to a good 
approximation) : 



eA ■ Pb ^ £B ■ PA ^ 0. 



(7) 



Finally, dot products of the form eA,B • P will be 
velocity-suppressed and can be safely neglected. 
Thus, of the ten original tensor structures which 
make up Al'"', only the 5'^" term survives. We 
have explicitly checked that, indeed, the Ai term 
is the dominant contribution even when all of the 
tensor structures are kept. 

In the calculation of the loop amplitudes, we ap- 
ply the following algorithm. First, we reduce all 
dot products of the form k ■ p (where k is the loop 
momentum and p is a generic external momentum) 
against the corresponding factors from propaga- 
tors in the denominator. The resulting three-point 
functions can be safely computed using the stan- 
dard Passarino-Veltman (P V) technique 2l| . The 
remaining terms in the amplitudes take the form of 
four-point scalar {Dq) and rank-two tensor {D^^) 
functions: 



Do- i_ii,{pi,P2,p3; mi, m2, 1713,7714) 
d''k 



{1; kfj^ki/} 



in^ _ ml][{k + pi)^ - ml][{k +pi +^2)^ - m|][(fc + p2 +^3)^ 



(8) 



where pi are external momenta and mi are the 
masses of the particles circling the loops. Note 
that we have neglected rank-one tensor integrals 
{Dfj_) since we are only interested in extracting the 
5^1^ pieces of the amplitude. 



Following the PV scheme, the rank-two tensor 
integral can be rewritten as a linear expansion in 
tensor structures which are built from the external 
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FIG. 2: Examples of Feynman diagrams which contribute to BrBh 'jV where V = -y, Z and B^^'^K 



momenta and the metric tensor: 

-D^i. = D2lPl,f_iP2.u + D22P2,iiP2.y H h -D275a'''' ' 

(9) 

The coefficients of this expansion (_Dy ) can then be 
reduced to scalar integrals [^|. However, in cases 
where two of the external momenta become iden- 
tical, as for the case of WIMP annihilation, this 
approach breaks down. In these cases, the expres- 
sions for the Dij coefficients in terms of scalar inte- 
grals depend inversely on the Gram Determinant 
(GD) built from the external momenta (i.e., GD 
= det{pi -pj)). In certain kinematical regions (e.g., 
where two of the momenta become degenerate and 
GD ~ 0), the PV scheme gives rise to spurious 
divergences. In calculations for collider processes 
(where the momenta are integrated over an entire 
phase space), this situation arises only at special 
points near the boundaries of phase space. Spe- 
cial techniques involving interpolating from these 
unsafe regions of phase space to safe regions have 
been developed to deal with these spurious diver- 
gences in calculations for collider processes. 

These techniques do not apply to our situation 
(where the two incoming momenta are fixed and 
identical) and we are forced to approach this prob- 
lem using the following method. For our calcula- 



tion, we have chosen to implement the technique 
developed in Ref. [12] ■ In this algebraic reduc- 
tion scheme, the original PV scheme is extended 
to deal with situations where the GD exactly van- 
ishes. Higher-point tensor (and scalar) integrals 
are expressed in terms of lower-point quantities 
which can be safely evaluated utilizing the usual 
numerical techniques. For example, the expres- 
sions for the four-point scalar integral (i^o) a-nd 
the tensor coefficient (-D27) can be expressed as: 



Do = ai23Co(123) 
+ ai34Co(134) 



ai24Co(124) 
a234Co(234) , (10) 



and: 



£'27 = ai23C24(123) + ai24C24(124) 

+ ai34C24(134) + a234C24(234), (11) 

where C'o{ijk) and C24:(ijk) are the three-point 
scalar integral and PV tensor coefficient, respec- 
tively (the (ijk) denotes various propagator factors 
in the original four-point denominator). The Uijk 
coefficients can be obtained by solving the matrix 
equation: 
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/II 1 1 

pI {pI-pI+pI)/2 {pI+pI-pI)/2 

{-pI-pI+pI)/2 [-pI+pI+pI)/2 [^pI-pI+pI+pI)/2 

2 9 2 2 2 2 2 

— mi PI — 7712 Pt ~ ^3 Pt ~ ^4 





/ "234 \ 






1 


^134 


- 






ai24 








\ q;i23 / 







(12) 



where pi , . . . , p4 are the external momenta, p^ = 
Pi + P2, P6 = P2 + P3 and the are loop particle 
masses. 

This approach allows us to construct quite 
compact expressions for the one-loop amplitudes 
for the processes of interest in terms of kine- 



matical factors and scalar integrals of the form 
Bo{p^;ml,ml) and Co{pl,pl,pl; ml, mj, ml). For 
example, the total amplitude for BhBh — > 77 via 
loops of one (massless) fermion species and its KK 
partners is given by: 



A 



-ayaemQUYl + Y^) 1^2 + ^Bo{Ml^ ; M|, 0) - i?o(4Af|^ ; 0, 0) - i±^i?o(4Af|^ ; Ml Ml) 
+Ml^ -{l + r^){Co{Ml^AMl^,Ml^;MlO,0) + Co{Ml^AMl^,Ml^-,0,MlMl)) 

l^; Ml Ml Ml) 



-2Co(A/|^, 0, Ml^;0, Ml Ml + AfiCo{0, 0, 4MI 



(13) 



where aem and ay denote the SM fine struc- 
ture constant and U{1) coupling constants, respec- 
tively. The charge of the fermion is (in units 
of e), while Il^/j are the left- and right-handed 
hypercharge quantum numbers of the fermion. In 
the above and the following, we assume a com- 
mon mass for all KK fermions Ml (i.e., we ne- 
glect any mass splittings) and we define rj = 
Ml/M^^. Similar (though more complicated) ex- 
pressions hold for the amplitudes for the processes 
BhBh iV where y is a massive gauge boson 
(e.g., a SM Z boson or a higher KK mode such as 
the B(i'i)). 

In Fig. [31 we plot the annihilation cross sections 
for 77, Z7 and B^^'^^^ production as a function 
of the LKP mass Mb„- Here and in what fol- 
lows, we have exhanged the compactification scale 
L for the LKP mass and then derived the other 
KK masses {Ml and Mg{i,i]) in terms of Mbh: 
using Ml = 1.17A/b„ and M^d.i) = 1.6Mb^, but 
our results are not strongly dependent on these as- 
sumptions unless the mass splittings are unusually 
large. It is interesting to note the enhancement 
of the B^^'^^j cross section compared to the 77 
and Z7 cross sections. This effect can be under- 
stood as follows. First, in the cases where V is 
either a photon or a Z boson, the couplings be- 



tween V££ and V^^ are nearly the same strength 
and the result is a significant cancellation between 
the various diagrams. However, in the case where 
V is identified with the S*^^'^) KK gauge boson, 
the Vi£ couplings are loop-suppressed, while the 
V^^ couplings are relatively large. This results in 
less cancellation between the various diagrams and 
an enhanced cross section compared to the other 
two processes. 



A. Line Spectra 

While the spectrum of the 77 line is simply a 
delta function at M Bh 1 the lines arising from the 
process Bh Bh — > 'jV will exhibit a line with an 
intrinsic width, which will depend on the mass of 
the boson in the final state My, 



U'Ib^MvT 



dE /1/2 
where Fy is the V width and 

Mv ■ 



(14) 



/: 



tan 



M 



+ tan 



4Af|^ - M^ 
MvTv 



f2 = 



M^y+T^yM^ 
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FIG. 3: Cross sections as a function of the LKP mass 
rnBu for the three different channels 77 (solid), -yZ 
(dotted) and -yB'-^''^'' (dashed). 



For Bh Bji Z7, we use the experimentally mea- 
sured values of Mz — 91 GeV and Tz — 2.5 
GeV m. For BhBh ^ B^^'^^-f, we have as- 
sumed the minimal mass boundary term relation 
Msd.i) « 1.6Mb„. The width of S^^'i) is deter- 
mined by its decays into SM particles. Such inter- 
actions result entirely from boundary terms, and 
under the minimal assumption that they are loop- 
suppressed [1^, we expect small coupling to the 
SM such that F^d,!) « lO-'^M^d,!) . Since F^d.i) 
is typically much smaller than the typical experi- 
mental energy resolutions, the resulting signals at 
detectors are not very sensitive to its precise value. 

Generally, the larger the mass of V , the further 
away its corresponding line will be from the lim- 
iting value Mbh- The chiral square model is in- 
teresting in this respect because the S'-'^^^) mass 
is large with respect to both Mz and , while 
still kinematically allowed. This leads to a pho- 
ton distribution for the Bh Bh B^^'^''^ process 
that peaks at energies that are far below the ones 
associated with the and 77 processes, which in 
turn makes the B^^'-^^^^ line clearly distinguishable 
from the other two lines even by current experi- 
ments, characterised by a relatively large energy 
resolution. 



V. PROSPECTS FOR DETECTION 

In order to predict the 7-ray flux from the galac- 
tic center generated by the chiral square model, 
we turn to the evaluation of the spectra for the 
Bh Bh 'yV processes and of the astrophysical 



uncertainties related to the integration over the 
dark matter distribution. The differential flux of 
photons arising from dark matter annihilation ob- 
served in a direction making an angle ■0 with the 
direction to the galactic center (GC) is given by 



dVLdE 



with 



dN^ 



4ttMI^ dE 



r ds 


>[r(s, V)]" 


/l.o.s. ^0 


Pq 




(15) 



dN-y 

He 



dNf 
dE 



(16) 



where we use the index / to denote the annihila- 
tion channels with one or more photons in the fi- 
nal state, {av) / is the corresponding cross section 
and dNj^ /dE is the (normalized) photon spectrum 
per annihilation. Furthermore, p{x), pq = 0.3 
GeV/cm'^ and tq — 8.5 kpc respectively denote 
the dark matter density at a generic location x 
with respect to the GC, its value at the solar sys- 
tem location and the distance of the Sun from the 
GC. Finally, the coordinate s runs along the line of 
sight, which in turn makes an angle ijj with respect 
to the direction of the GC. 

The specification of the dark matter profile is 
the largest source of uncertainty in the evaluation 
of the detectability of a dark matter annihilation 
signal arising from the galactic center, as it fixes 
the normalization of the predicted flux. The most 
recent high-resolution numerical simulations show 
that DM halos can be reasonably well fit with the 
Navarro Frenk and White (NFW) profile, which 
is often used as a benchmark for indirect searches 

M 



(17) 



However, modifications of the above profile on very 
small scales have been observed in the most re- 
cent simulations. In Ref. [25( it was argued that 
the innermost regions of DM halos are better ap- 
proximated with r~^'^ cusps, while in Ref. [26^1 it 
was found that the analytic form that provides an 
optimal fit to the simulated halos is the so-called 
"Einasto profile" ^ 



p{r) = po exp 



(18) 



which is shallower than NFW at very small radii^ . 



^ The values assumed for the parameters are in this case 
a = 0.17 and _R = 20 kpc. 
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Model 


J{W-'') 


NFW 


1.5 X 10" 


Adiabatic 


4.7 X lO'^ 



tered on the GC 



TABLE I: Value of J(10"^) for the dark matter density 
profiles adopted in Fig. |4] 



These results have been derived in the frame- 
work of simulations containing only dark matter, 
not taking baryons into account, which are ex- 
pected to play an important role in the dynam- 
ics of galaxies, especially on small scales. Due to 
dissipative processes, baryons in fact lose energy 
and contract, thus affecting the gravitational po- 
tential experienced by DM. In the "adiabatic com- 
pression" scenario , the baryons contraction is 
quasi-stationary and spherically symmetric. If one 
starts from an initial NFW profile, the final slope 
in the innermost regions is r~^'^ 0, H^, IMI) i32 |. 
Other processes, such as the "heating" of the DM 
fluid due to gravitational interaction of baryons, 
may have the opposite effect on the final DM dis- 
tribution, thus depleting the central cusps. 

Furthermore, the presence of a super-massive 
black hole (BH) will inevitably affect the DM pro- 
file, especially within its gravitational radius (i.e. 
where the gravitational potential is dominated by 
the black hole itself). The growth of the BH from 
an initial small seed would initially lead to a large 
DM overdensity, or 'spike' ^3§\. Dynamical effects, 
and DM annihilations, will subsequently tend to 
destroy the spike [H, [H, [11], although in some 
cases significant overdensities can survive over a 
Hubble time. In what follows, we will show our 
results for two profiles: NFW, and the adiabati- 
cally contracted profile, with the same parameters 
as the profile labelled 'A' in Ref. 

In order to quantify the astrophysical uncertain- 
ties related to the dark matter density profile and 
to separate them from the microphysics, we iden- 
tify with J the dimensionless integral along the 
line-of-sight appearing in Eq. and with J its 
average value computed for a solid angle Afi cen- 



^ Note that our values of J for the NFW and "Adiabatic" 
profiles are slightly different from the values in Ref. [33 , 
due to the fact that the NFW profile was there approxi- 
mated as a simple power-law from the galactic center 
to the location of the Sun. 



r ds 


'p[r{s,ip)]' 


h.o.s. ''0 





JiAn) = — 



j{'4j)dn. 
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In Tab. [J we show the values of J for the two 
halo profiles adopted in Fig. 21 obtained assum- 
ing Afi — 10^^ sr., corresponding to the angular 
resolution of the HESS and Fermi LAT 7-ray ex- 
periments. 

To account for the finite resolution of the detec- 
tor we convolve the unfiltered "raw" signal S{E) 
with a gaussian kernel G{E, Eq), 



GiE,E„) 



1 



2'kEq(t 



exp 



{E~Eo 



i21 



2(7^ El 



, (20) 



where a is related to the detector's relative energy 
resolution ^ by ct = 'C/2-3. The signal Sm{Eq) 
measured by the detector at energy Eq is then sim- 
ply given by 



S 



(Ea) - J dEG{E,Eo)S{E). 



(21) 



We show in Fig.[3]the final predicted fluxes, from 
a solid angle Afi = 10~^ towards the galactic cen- 
ter, for the chiral square model with Mbh = 250 
GeV and Mb^ = 450 GeV and for two different 
choices of the DM profile. Masses of relevant states 
for these two mass scenarios are given in Tab. [Hi 
We show both the actual spectrum (dotted lines) 
and the spectrum as it would be observed by an ex- 
periment with a 10% energy resolution (solid) such 
as Fermi-LAT. We also show for comparison the 
HESS data for the point source HESS J1745-290, 
which is spatially coincident, within the angular 
resolution of HESS, with the supermassive black 
hole Sgr A* lying at the center of the Galaxy [s^ . 



VI. DISCUSSION AND CONCLUSIONS 

As can be easily seen from Fig. HI the WIMP 
forest in the chiral square scenario is characterized 
by three lines. The two lines at and around Mb„, 



Mb„ 




Ml 




250 


290 


294 


401 


450 


521 


528 


721 



TABLE II: Masses of Bh , Ci^l and B'^'^' and curvature 
in GeV for the two mass scenarios given in Fig. |4l 
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FIG. 4: Predicted fluxes, from a solid angle AQ — 10~^ towards the GC, for the chiral square model with 
Mbh = 250 GeV (left column) and Mbh = 450 GeV (right). We show both the actual spectrum (dotted lines) 
and the spectrum as it would be observed by an experiment with a 10% energy resolution (solid) like Fermi 
LAT. An NFW (adiabatically compressed) profile has been adopted for the lower (upper) panels. We show for 
reference the HESS data relative to the gamma-ray source detected at the Galactic center. 



belonging respectively to the Bh Bh —>■ 77 and 
the Bh Bh Z-f processes cannot be resolved by 
a detector whose resolution is only 10% and they 
appear as a single bump. On the other hand, the 
Bh Bh B'-^'^'>^ process is really the distinctive 
feature of the chiral square model, leading to a line 
at much lower energies that is clearly distinguish- 
able from the other ones even by detector with only 
fair energy resolution. The information content of 
such a line can be hardly overemphasized. Since 
5D Kaluza-Klein models are lacking such a fea- 
ture, experimental observation of this line would 
not only make a strong case for Kaluza-Klein dark 
matter but also indicate the presence of two extra 
dimensions compactified along the lines of the chi- 
ral square model and allow a first measurement of 
the size of such extra dimensions. 

Figure m also shows that, for a suitable choice of 
astrophysical and particle physics parameters (es- 
pecially for an 'optimistic' DM profile), the anni- 



hilation flux in this scenario is within the reach 
of current experiments. The size of the error 
bars in the HESS data, also suggests that the 
WIMP forest (an abbreviated name for the "forest 
of gamma-ray lines produced by WIMP annihila- 
tions" ) might be realistically detected as a "double 
bump" in the energy spectrum. It should be kept 
in mind, however, that the HESS source actually 
exhibits a flat, feature-less, spectrum extending 
over two decades in energy. This means that the 
underlying emission probably originates from an 
ordinary astrophysical source, and that this source 
represents a foreground for a DM-induced signal. 

A better strategy might consist in the analy- 
sis of the gamma-ray emission from a broader re- 
gion, e.g. a cone of several degrees towards the 
center, or an annulus that excludes the galactic 
plane [s^, [s^. A dedicated analysis of the var- 
ious detection strategies, also containing a dis- 
cussion of the prospects for detecting gamma-ray 
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lines with the Fermi LAT, has been recently pub- 
lished in Ref. [37[, where the authors conclude 
that for robust identification, annihilation cross- 
sections to lines as high as av ~ 10~^^ cm^ s~^ are 
required, for an NFW profile. With steeper pro- 
files, even smaller cross sections lead to annihila- 
tion fluxes comparable with, or larger than, astro- 
physical backgrounds, thus becoming potentially 
accessible to gamma-ray telescopes. In Ref. ^SQi] , 
for instance, it is shown that a suitable choice of 
the observational window may allow to probe cross 
sections of order av ~ 10~^^ cm^^ s~^ for profiles 
steeper than r~^'^ . These results however strongly 
depend on the intensity and shape of astrophysical 
foregrounds, that the authors approximated with 
an analytic fit to EGRET data proposed over a 
decade ago in Ref. [i^. Preliminary Fermi mea- 
surements suggest that EGRET may have overes- 
timated the flux of photons at GeV energies, there- 
fore a word of caution is in order when estimating 
the prospects for detecting gamma-ray lines in this 
energy range. 

Several other targets for indirect detection have 
been proposed over the years (see Refs. [l| for 
a discussion and references). The current under- 
standing is that, aside from a small region of size 
C(l) degrees from the galactic center, the diffuse 
flux is dominated by annihilations in Galactic un- 
resolved substructures [4l|, [13], while the extra- 
galactic contribution is expected to be subdomi- 
nant [11]. Alternative targets include clumps of 
DM (e.g. [13, Si], H^l and references therein) and 
mini-spikes around Intermediate Mass Black Holes 



In closing, gamma rays from the annihilation of 
dark matter are a fascinating signal that can tell 
us a lot about the underlying particle theory, and 
probe the distribution of WIMPs in the galaxy. 
The lines of the WIMP forest serve to measure the 
mass of the WIMP, detect the presence of other 
states in the dark sector, and perhaps even iden- 
tify which theory of dark matter is responsible for 
a signal. While we have worked in the specific con- 
text of a six dimensional UED model, the feature 
itself is more general, and can be present any time 
two WIMPs can annihilate into a photon and an- 
other massive boson. As we enter an exciting new 
era of possibilities to detect WIMPs indirectly, it 
is important to search for multiple features in the 
spectrum of gamma rays. 
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